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ABSTRACT 



CD 

■ Distribution of galaxies may be a biased tracer of the dark matter distribution and 

the relation between the galaxies and the total mass may be stochastic, non-linear 
. and time-dependent. Since many observations of galaxy clustering will be done at 

^!t| ! high redshift, the time evolution of non-linear stochastic biasing would play a crucial 

role for the data analysis of the future sky surveys. Recently, analytic study of the 



(N 



Oh' 



0^ I time evolution induced by gravity has been reported in the mildly non-linear regime. 

' Here, we further develop the non-linear analysis including the next-to-leading order 

On ■ and attempt to clarify the non-linear feature of the stochastic biasing. Employing the 

perturbative approach, we compute the one-loop correction of the power spectrum 
for the total mass, galaxies and their cross correlation. Assuming that the initial 
2 , distribution of galaxies is given by the local function, we specifically investigate the 

I time evolution of the biasing parameter and the correlation coefficient deduced from 

the power spectra. On large scales, we find that the time evolution of the biasing 
parameter could deviate from the linear theory prediction in presence of the initial 
^ ' skewness, even though the scale-dependence of the biasing is very weak. On the other 

■ hand, the deviation can be reduced if the stochasticity between the galaxies and the 

total mass exists. To explore the influence of non-linear gravity, we focus on the 
quasi-linear scales, where the non- linear growth of the total mass becomes important. 
It is recognized that the scale-dependence of the biasing dynamically appears and the 
initial stochasticity could affect the time evolution of the scale-dependence. The result 
is compared with the recent N-body simulation that the scale-dependence of the halo 
biasing can appear on relatively large scales and the biasing parameter takes the lower 
value on smaller scales. Qualitatively, our weakly non-linear results can explain this 
trend if the halo-mass biasing relation has the large scatter at high redshift. 



Subject headings: cosmology:theory— large scale structure of universe, galaxies: biasing 
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1. Introduction 

Understanding the large scale structure of our universe is the ultimate goal of observational 
cosmology. Currently, Sloan Digital Sky Survey and Two-Degree Field Survey are working and the 
enormous observational data will be obtained. In the next generation of deep galaxy surveys, the 
formation and the evolution of galaxy distribution will be measured more precisely. It provides us 
with the powerful constraint on the theory of structure formation in the Universe. 

When we compare the theoretical prediction for the density fluctuations with the observation 
of galaxy sky survey, however, the relation between the total mass and the galaxies must be 
clarified. Since we do not yet have any reliable theory of galaxy/star formation history, there exists 
the statistical uncertainty between the galaxies and the total mass, which must be determined 
from the observational data. This uncertainty hampers the effort to obtain the cosmological 



information such as the density parameter (Hamilton 1998). 



Provided the fluctuations of the total mass 6m and the distribution of the galaxy number 
density 5g, the relation 

5g = b 5m (1) 

is often quoted in the literature. However, the galaxy biasing in nature could be non-linear and 
the simple linear relation (|^) might not be appropriate ( Fry &: Gaztahaga 1993| ). Fry (1994) and 



Gaztahaga & Frieman (1994) compared the observational data with a straightforward extension 
of (|l|) given by 

= f{5m) = E ^ i^rnT- (2) 
n=l 

The assumption (pi) still restricts the statistical feature of the galaxy distribution. In general, the 



relation between the galaxies and the total mass may be time dependent and stochastic (Cen & 



Ostriker 1992| ). Therefore, a more general framework to treat this situation is needed, referred to 



as the stochastic biasing. 

Using the different catalogs or subsamples of the galaxy surveys, the observational evidence 
of stochasticity has been recently discussed ( [Tegmark &: Bromley 199S , Seaborne et al. 1999). As 



for the theoretical study, Dekel & Lahav (1999) has proposed the general formalism for non-linear 
stochastic biasing. Pen (1998) discussed an observational method to determine the stochasticity. 
Further, Blanton et al. (1998) and Blanton, Cen & Tegmark (1999) performed the cosmological 
simulation of galaxy formation and demonstrated that the stochastic property of the biasing could 
evolve due to the galaxy formation process. Since we still lack our knowledge of the non-linear 
aspect of the stochastic biasing, the qualitative features should deserve consideration . 

In this paper, we analytically investigate the dynamical feature of the non-linear stochastic 
biasing. We study the time evolution induced by gravity, which would be one of the most 
important sources of non-linearity. The observational study of the galaxy clustering recently 
shows that the evolution model induced by gravity is sufficient to explain the observed galaxy 
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distribution at redshift z < 2 (Magliocchetti et al. 1999, Magliocchetti et al. 199(^ ). Therefore, the 



gravitational evolution of biasing would play a crucial role in the data analyses of future deep sky 
surveys. 

The previous discussion about the time evolution of stochastic biasing has been concentrated 
on the tree-level analysis, i.e, based on the lowest order perturbation. Tegmark & Peebles (1999) 
examined the linear evolution. Taruya, Koyama &: Soda (1999) has extended the pioneering work 
of Fry (1996) and investigated the skewness and the bi-spectrum. Taruya & Soda (1999) further 
study the weakly non-linear galaxy-mass density relation. 

This paper focuses on the next-to-leading order of non-linearity. We take into account the 
one-loop correction of non-linear gravity, corresponding to the one-loop diagram in the graphical 
representation of perturbation theory ( pcoccimarro &: Frieman 1996a , 1996b). Then, we shall 



investigate the power spectrum, the second order statistics of density perturbations in the Fourier 
space. The power spectrum of the total mass distribution has been extensively investigated by 
several authors ( [Makino, Sasaki fc Suto 1992| , Jain &; Bertschinger 1994, Baugh & Efstathiou 



1994, Scoccimarro &: Frieman 1996b| ). Our analysis is further extended to the power spectrum of 



galaxies and the cross correlation between 5m and 6g. Define the vector X'^ = {6m{k), 5g{k)) in the 
Fourier space, we have the three kinds of power spectrum 

= (2.f^o(>^. + ( ^;<,^;> ) . (3) 

where Pm{k), Pg{k) and Px{k) are the power spectrum of the total mass, galaxies and the cross 
correlation between them. The statistical feature of the galaxy biasing is encoded in the power 
spectra Pg{k) and Px{k), which can be quantified by comparing them with the power spectrum 
of the total mass distribution. We will analyze the biasing parameter and the correlation 
coefficient : 

h - r m 

V ^'»(^) ^Pg{k)Pm{k) 



which are frequently used in the literature ( p?egmark fc Peebles 199^ , [Tegmark &: Bromley 1999 ). 



Note that the simple biasing prescription (|l]) restricts the above parameter r/j = 1. 

We organize the paper as follows; In section |2|, we briefly mention the basic formalism for 
the time evolution of stochastic biasing. The evolution equations (eqs.[^0 |§) and the initial 
conditions (eqs.fl^ |jl^) are presented. Section ^ is devoted to the weakly non-linear analysis. We 
compute the one-loop correction of the power spectra and the time evolution of the stochastic 
biasing is investigated in detail. Our main results are described in Section We specifically pay 
attention to the effect of the non-linear biasing (Sec. 3.3.T| ) and the non-linear gravity (See p. 3.2 ). 



The conclusions and discussion are summarized in section 
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2. Gravitational evolution of Stochastic biasing 

Recall the definition of the density field 5m and 5g: 

Smix) ^ 5,i.) = !^MZ^, (5) 

where the variables with overbars denote the homogeneous averaged quantity. 

The gravitational evolution of the galaxy biasing has been studied by Fry (1996), Tegmark & 
Peebles (1998) and Taruya, Koyama Sz Soda (1999). In these papers, the distribution 6m follows 
the equation of continuity and the Euler equation. On the other hand, neglecting the merging 
and the galaxy formation process, the galaxies are regarded as the test particles in the dark 
matter distribution Sm- Thus, the distribution 5g should satisfy the equation of continuity whose 
velocity field is determined by the gravitational potential. On large scales, the assumption of the 
irrotational velocity flow is valid. We define the velocity divergence = V • v/{aH), where a is the 
scale factor of the universe and H is the Hubble parameter. The evolution equations for the total 
mass and the galaxies become 

^ + He + -V-{6mv) = 0, (6) 
at a 

f + «« + iv^(M = o. (8) 

The variable A is the cosmological constant and 17 is the density parameter defined by 

n = ^ A (9) 

Equations (0) and ^ are our basic equations for the time evolution of 6m and 6g. 
These equations can be solved systematically using the perturbation theory if we give the initial 
conditions 6m,g{x,ti) at an early time tj order by order. To linear order, provided the initial 
distributions 6m and 6g are given in terms of the independent random fields /S.m{x) and /S.g{x), the 
solutions become 

5«(x) = i?(t) A„(x), 6^^\x) = {D{t)-l}^m{x) + l^g{x). (10) 

Since we are interested in the leading behavior of the time evolution, the solutions are only taken 
into account the growing mode of perturbation denoted as D{t). The function D{t) satisfies 

ij + 2HD -^H^nD = 0, (11) 

with the initial condition D{ti) = 1. Note that the random fields and assign the statistical 
feature of the galaxies and the total mass. On large scales, the linear perturbation for the 
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fluctuations 6m]g may be valid and the central limit theorem implies that the spatial distribution 
could be well-approximated by the Gaussian statistics. Hence, the variables and are 
assumed to be the Gaussian random fields. 

The full evolution is non-linear, including the higher order perturbative solutions. Regarding 
the spatial variables A^.g as the seeds of perturbation, the general form of the n-th order solutions 
becomes the terms collecting the time dependent function multiplied by the spatial random 
variables [Am,g]". Here, we give the initial distributions 6m, g as the Taylor expansion of the 
functions (|Taruya, Koyama Sz Soda 1999| , Taruya &: Soda 1999 ) 



Sm{x,ti) = /(A„), 6g{x,ti) = g{Ag). (12) 



Based on the ansatz (12), we will consider the time evolution emanating from the local biasing. 
Their functional forms are determined as follows; since the total mass fluctuation is produced 
during the very early stage of the universe, the gravitational instability induces the deviation from 
the initial Gaussian distribution. Thus, we give the initial condition of the total mass /(A^) 
iteratively from the perturbative solution by dropping the decaying mode. As for the galaxies, the 
initial distribution 5g is induced by the galaxy formation. Currently, it is not feasible to know the 
initial non-Gaussian distribution g(Ag). We treat g as a parametrized function, 

<7(A,) = A, + ^(A^-(A^)) + |a3 + .... (13) 

Neglecting the higher order terms, ( [T^ ) gives the same initial condition as in the linear order 
solution 

Notice that the resulting expression for the total mass SfYi leads to the non-local and 
non-Gaussian distribution (see eqs.[ [A-l| [ |A-3| ). On the other hand, we simply assume that the 
function g{Ag) is locally characterized by the constant parameters hi and /12, which have to 
do with the initial skewness and initial kurtosis. We should remark that the locality of the 
galaxy distribution only holds at the initial time. Due to the gravity, the distribution 6g can 



become non-local, which can yield the scale-dependence of the galaxy biasing (Sec. 3.3.2 , see also 



Matsubara 1999D 



Once we obtain the perturbative solutions 5m,g, we can investigate the statistical features for 
the density fluctuations by taking the ensemble average (•••). In our prescription, the ensemble 
average is taken with respect to the Gaussian random fields Am,g. Their stochasticity is completely 
characterized by the three quantities. Defining the two-dimensional vector = (Am, A^) in the 
Fourier space, we have 

(y(fcl)y\fe2)) = {27Tf6D{kl + k2) Plih) ( ^ ^ ] , (14) 



where the quantity Pi{k) is the initial power spectrum of the total mass. The parameters 60 and 
ro are assumed to be the scale-independent variables, for simplicity. 
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Now, the formalism developed here enables us to obtain the weakly non-linear power spectra 
@. Before closing this section, we present the linear order results. Using the relation (p!4[), the 
three kinds of the power spectra are obtained by substituting the Fourier transform of the linear 
solutions ( p!o[ ) into the definition (|3|). Denoting the linear spectrum as P^^^\k), we have 

Pjj''\k) = {{D - if + 2boro{D - 1) + bl} Pi{k), (15) 
P^^'\k) = D{iD-l)+boro} Pi{k). 

Thus, the biasing parameter bk and the correlation coefficient can be evaluated from ^ 
( Tegmark fc Peebles 199S , Taruya, Koyama &: Soda 19991) : 



where the superscript ((')) means the linear order quantity. The variables 6^'^ and r^') are 
scale-independent and their initial values correspond to the parameters 6o and tq. Note that Fry 
(1996) has obtained the same result in the specific case ro = 1. In this case, we immediately obtain 
r(^) = 1 and the relation between Sm and 5g reduces to the deterministic biasing (|l]). Thus, in our 
prescription of the time evolution, the stochasticity in the galaxy biasing comes from the initial 
condition tq < 1. 



3. Perturbative approach to the non-hnear power spectra 

In this section, we develop the higher order perturbations and investigate the time evolution 
of galaxy biasing. In section 3A, the one-loop correction of power spectra is computed. After 
describing the qualitative feature of the galaxy biasing in the weakly non-linear regime, the time 
evolution is examined in section 3.3. 



3.1. One-loop corrections of power spectra 

When we take into account the next-to-leading order terms, the power spectrum can be 
written as 



Pim,,,.){k) = PLl.)(k) + ^^^Ix) W + ^(^Ix) W + ■ ■ ■ • (17) 



5(22) ^ p(13) 



The terms in the bracket are the one- loop correction, which can be evaluated using the second 
and third order solutions of perturbations Sg and 5m- The details of calculation and the resulting 
expressions for the one-loop spectra P^'^'^\k) and P^'^'^\k) are comprehensively described in 
appendix. Here, we discuss the technical issue to implement the numerical integration. 
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The expressions for the spectra ( A-7[ )-( |A-12 ) have the three-dimensional integral for the 



linear power spectra Pi{k). The further reduction enables us to obtain the expressions including 
the two-dimensional integral over u and x for the P^'^'^\k) part, while the P^^'^\k) part has 
the one-dimensional integral over x (see appendix). The care is required for evaluation of these 
integrals. For the single power-law spectrum Pi{k) oc fe" with n < —1, the infrared divergence 
appears at each integral. As for the total mass Pm (k) and Pm {k), which have been studied by 
several authors, the limit x — > yields 

Pj^^\k), P^^\k) OC ePi{k) / dqPi{q), 

Jo 

which imply that the infrared divergences occur at n < — 1. It is known that this is apparent 
and the divergences cancel out when we evaluate the total contribution of the power spectrum 
Pmik) ( [Makino, Sasaki fc Suto 1992| , |Jain Sz Bertschinger 1994 ^coccimarro &i Frieman 1996bD. 



However, in cases of the power spectra Pg{k) and Px{k), the infrared divergences still remain even 
if we compute the total contribution. 

The physical interpretation of the infrared divergence in cases with the single power-law 
spectrum Pi{k) is difficult. In reality, however, the infrared divergence does not appear if we adopt 
the phenomenological spectrum such as the Cold Dark Matter (CDM) power spectrum by Bardeen 
et al. (1986), where the effective spectral index ?T-e// = dlogPi{k)/dlogk becomes positive in the 
limit /c — > 0. Therefore, the particular behavior in the single power-law spectrum might not affect 
the real universe. 

On the other hand, depending on the shape of the power spectrum Pi{k), the ultraviolet 
divergence (x oo) appears even in the total mass Pm{k) ( [Scoccimarro fc Frieman 1996b|). The 



integrand from the high frequency part can give the large contribution to the total power spectra 
and the care is also required to proceed the weakly non-linear analysis self-consistently. Thus, 
the high-frequency cutoff kc should be introduced when we evaluate the integrals. If we have the 
power spectrum Pi{k) as defined with the interval e < k < kc, the integral over x and u for the 
P(22)(/^) part and over X for the P^^^\k) part should be replaced: 



oo 



kc/k rl |-Min{l,(l+X'2-(fce/fc)2)/(2x)} 

dx — > / dx, / du — > / du. (18) 



Je/k J~l JMa.yi{-l,{l+x^-{e/k)'^)/{2x)} 

The cutoff e is auxiliary introduced so that the numerical integration can converge rapidly. As long 
as the power spectrum Pi{k) has the effective index n^jf > —1 in the limit A; ^ 0, the influence of 
the cutoff e does not affect the weakly non-linear power spectra. 



3.2. Non-linear biasing and power spectrum 



The qualitative features of the one-loop power spectra should be mentioned before analyzing 
the time evolution. Several authors have discussed the influence of non-linear biasing to the 
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two-point correlation function ( 


Coles 1993, 


Fry & Gaztanaga 1993, Scherrer & Weinberg 1998, 


Heavens, Matarrese & Verde 1998 




Dekel & Lahav 1999 


). The remarkable features appearing in 



the spectrum of galaxies are summarized as follows : 

(i) Constant power on very large scales, 

(ii) Over a wide range of scales, the non-vanishing parameters hi and /12 modify the biasing 

parameter predicted in linear theory. 

The effect (i) comes from the Pg {k) part of the power spectra. From ([A- 91), the initial 



one-loop spectrum Pg 
hi 



^'^'^\k) becomes 

2,3 

1 tXi 



18 " (27r)2 



For the typical linear spectrum Pi{k) satisfying the limit Pi{k) — > as /c ^ 0, the above integral 
becomes constant and it eventually dominates over both the linear term Pg^^^ and Pg^^^ part. 
This leads to the divergence of the biasing parameter on very large scales. 

(22) 

According to Heavens, Matarrese & Verde (1998), the constant power in the Pg part is 
interpreted as the shot-noise term, which is a consequence of the local ansatz for the galaxy 
distribution (see eq.||l^). Although the constant power can be measured on extremely large scales, 
this would become negligible on realistic scales of the galaxy surveys. Therefore, we will analyze 
the power spectrum Pg{k) on the scales where the effect (i) becomes negligible. 

On the other hand, the clustering of galaxies incorporated by the non-Gaussian initial 
condition ( [T^ leads to the effect (ii), which would have a more practical importance for the data 
analysis of galaxy surveys. The dominant contribution of this effect comes from the P^^^\k) part. 
In our case, the initial biasing parameter is estimated using equation ( [A-10| ): 



\ 



pr^ + p 



(13) 



P, 



(11) 



(2^ 



(19) 



The expression (|T^ shows that the deviation from linear theory prediction depends on the 
parameters do, 60 and /12. The spectrum Pi{k) with the cutoff parameter implies that the 
parameter ctq would have the cutoff dependence. In a realistic situation, the clustering property 



characterized by the assumption (13) might in general depend on the scale. Hence, it is natural to 
have the cutoff dependence. As for the initial correlation coefficient, it remains unchanged due to 
the assumption (^). The expressions (A-10) and (A-12) yield 



rk 



,(11) 

X 



,(13) 

X 



(pr + p^); 



Pr 



(11) 



1 _L hz h2 2 
i -I- g OqCTq 



1 _L M fj2 2 



ro + Oial), 



which is valid as long as fio ^ 1. 
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Notice that the gravitational evolution alters the effective biasing and the correlation 
coefficient. The non- linear biasing parameters hi and /12 of the one-loop spectra dynamically 
affect these parameters. Further, the non-linear growth of density fluctuations might change the 
evolution. 



3.3. Time evolution of and 

We now consider the time evolution of the biasing parameter and correlation coefficient 
Tfc. In the rest of this paper, we present the results for the linear spectrum Pi{k) given by the 
standard CDM model. Although we have also examined the two-power-law model of Makino, 
Sasaki & Suto (1992), we obtained the qualitatively similar results to those shown below. Thus, 
our results would be general in the case of the realistic spectrum with a break. 

From Bardeen et al. (1986), the fitting form of the linear spectrum Pi{k) becomes 

Pi{k) = Ak T\k), 
C 

2Mq 



T{k) = ^"^Vo!"^^'^^ + 3-89g + (16.1g)2 + (bA&qf + {Q.llqf ^'^ 



k 



r h Mpc" 



-1 ' 



(20) 



where k is in unit of h Mpc~^. We have chosen the shape parameter F = Q,Qh = 0.5, 
corresponding to the parameters VLq = 1 and Hq = lOO/i km s~^Mpc~^ with h = 0.5. The linear 
growing mode D{t) is therefore replaced with a{t)/a{ti). The cutoff parameters are specified as 
e = 0.0001 h Mpc~i and kc = 27r/8 h Mpc-\ appropriate for examining the weakly non-linear 
evolution on large scales. Then, the amplitude A is normalized by erg = 0.54 at the present time 
(z = 0), determined from the cluster abundance ( [Kitayama fc Suto 1997] ). We then perform the 
numerical integration for the one- loop corrected terms. 



The evolved results of the weakly non-linear power spectra (17) are shown in Figure 1. 
The thin-solid, thick-solid, and thick-dashed lines are the power spectra Pm{k), Pg{k) and 
Px (k) evaluated at present, respectively. The initial conditions are set at z = 5. Choosing the 
non-Gaussian parameters hi and /12 as zero, we specify the parameters bo = 3.0, ro = 0.8 so that 
the initial biasing parameter and initial correlation coefficient obtained from (^) become bk = 3.0, 
rk = 0.8. 

In Figure 1, the evolved spectra for galaxies and cross correlation become the larger power 
than that of the total mass due to the initial large biasing. We have also displayed the linear 
power spectrum of the total mass Pm^\k) (see eq.||l5|). Compared to the linear order result, 
the weakly non-linear spectrum shows the enhancement of the power on small scales, while both 
Pm^\k) and Pm{k) have the same power on large scales. This behavior is known as the weakly 



non- linear effect ( Makino, Sasaki fc Suto 1992 , Jain &: Bertschinger 1994 ). The validity of the 
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weakly non-linear result has been checked using the non-linear fitting formula by Peacock & Dodds 
(1996). The one-loop correction of the total mass rather overestimates the non-linear growth at 
z = 0. However, until z = 1, the weakly non-linear spectrum Pm{k) is well-fitted to the formula 
by Peacock & Dodds and the qualitative behaviors did not change at z = 0. 

In Figure 2, using the definition (|^, we evaluate the biasing parameter bk and the correlation 
coefficient as the function of the Fourier mode k. Within the validity of the perturbation, the 
snapshots of both parameters with the same initial conditions as used in Figure 1 are depicted 
for the different redshift parameters z = 5,3, 2 and 1. On large scales, the spatial dependence of 
the biasing parameter and the correlation coefficient is very weak and they approach the constant 
values. The same behaviors can be obtained in the presence of initial non-Gaussianity, consistent 
with the results by several authors (|Mann, Peacock fc Heavens 199^ , pcherrer Sz Weinberg 1998 , 
Heavens, Matarrese Sz Verde 199^ ). On the other hand, the scale-dependence of the galaxy biasing 
exists on small scales. Although bk and initially lose their power due to the high frequency 
cutoff kc, the spatial variation of the biasing parameter and the correlation coefficient gradually 
changes. 

In the following, we shall separately analyze the time evolution in the linear scale and the 
quasi-linear scale. Here, the linear scale means that the evolution of the total mass power spectrum 
can be well- approximated by the linear perturbation (|l5|). The quasi-linear scale denotes the scales 
where the growth of non-linear correction for the total mass spectrum can become significant, 
corresponding to A; ^ 0.1 h Mpc~^ in Figure 1. 



3.3.1. Linear scale 

Since the non-locality of the gravitational evolution does not affect the biasing feature in the 
linear regime, the scale-independence of b^ and simply relies on the assumption of the local 



initial conditions (see eq.[13|). We thus analyze the time evolution by fixing the mode k = 27r/100 h 
Mpc~^. For our purpose, it is important to clarify the influence of the initial non-Gaussianity 
and the stochasticity to the evolution of bk and r^. Hence, by varying the non-Gaussian initial 
parameters hi and h2, we consider the two cases : almost deterministic biasing {bk = 3.0 = 0.8 
at z = 5) and stochastic biasing {bk = 3.0 Vk = 0.2 at z = 5). Because the initial non-Gaussianity 
slightly affects the initial conditions bk and r^, we appropriately adjust the model parameters b^ 
and ro so that the above conditions are satisfied. The list of the model parameters examined in 
this subsection is shown in Table 1, together with the evolved results bk{z = 0) and rk{z = 0). 
They are compared with the linear prediction ( [T^ ) under the same initial conditions bk and r^. 
The resulting evolution bk and is depicted in Figure 3a for the nearly deterministic biasing and 
Figure 3b for the stochastic biasing case. The results are naively extrapolated to the present time 
z = 0. 



First, consider the time evolution of the biasing parameter bk. In the deterministic biasing 
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case (upper panel of Figure 3a), the solid line shows the result with the vanishing initial 
non-Gaussianity. On the scale k = 27r/100 h Mpc~^, there still exists the scale-dependence of 
the biasing, as shown in Figure 2. This slightly alters the time evolution of bk, compared to the 
linear prediction 6^'^ (thin solid line). However, this deviation is small and the weakly non-linear 
evolution is in good accordance with the linear prediction until z = 1. The same behavior is 
obtained in the non-Gaussian initial condition /12 = 3.0 (dotted line). Although the non-Gaussian 
parameter /12 modifies the initial biasing parameter bo = 2.64 predicted by the linear theory (see 
eq. ) , the resulting evolution is well approximated by the linear evolution 6^') with the observed 
initial condition b^ = 3.0 and = 0.8. 

However, the deviation from the linear evolution is found in the presence of the initial 
non-Gaussianity hi ^ 0. The long-dashed and the short-dashed lines are the results with the 
initial parameter hi = 3.0 and hi = —3.0, respectively. In the one-loop power spectrum Pg (k), 



the terms proportional to the parameter hi has the time dependence D {t) (see eq.[A-10|), which 
is the same growth rate as appeared in the linear spectrum -Pj^^'* (see eq.[|l5|). This means that 
the clustering property induced by the initial non-Gaussianity tends to remain and the deviation 
from the linear theory would become increasing. In the /ii > case, we have the positively skewed 
distribution of galaxies, while the total mass distribution is nearly Gaussian. Thus, the fluctuation 
5g correlated with the total mass becomes more and more concentrated on the total mass, which 
dynamically leads to the larger biasing parameter bk. Vice versa, the smaller value bk could be 
attained in the case of the initial skewness hi < 0. 

The results of the stochastic biasing case {b^ = 3.0, = 0.2 at z = 5) are displayed in 
Figure 3b. For the biasing parameter b^ with the initial non-Gaussianity hi 7^ 0, remarkably, the 
deviation from the linear theory prediction is reduced. From the expressions (A-9) and ( |A-10| ), 



one can show that the one-loop corrected terms proportional to the parameter hi depends on 
the multiplicative factor 60^0 • The small value of the initial correlation coefficient — tq forces 
the galaxy distribution to correlate with the total mass. This behavior also affects the galaxy 
clustering and it can weaken the effect of initial non-Gaussianity. Hence, the biasing parameter bk 
could be approximated by the linear theory 6^'-*. 

Next, we consider the correlation coefficient. In the bottom panel of Figure 3, we plot time 
evolution of with the same initial parameters as examined in the upper panel of Figure 3. As 
was remarked by Dekel Sz Lahav (1999), the non-linear biasing such as the initial condition (|l^ ) 
can bend the 6g-5m relation and it effectively decreases the correlation coefficient, compared to 
the linear result In both panels of Figure 3a and Figure 3b, the effect of non-linear biasing, 
corresponding to the non-vanishing parameters hi and /12 is very weak. The time evolution of 
can be well-approximated by the linear prediction with the initial conditions bk = 3.0, = 0.8 or 
6fc = 3.0,rfc = 0.2. 



We have also analyzed the time evolution for more various initial conditions. The examples 
for the initial anti-biasing case are listed in Table 1. We see that the influence of the initial 
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non-Gaussianity to the parameters bk and becomes weak and the hnear prediction does work 
weh. In Table 1, the cutoff dependence kc is also examined. The small value of the cutoff 
wavelength 2-K/kc = 4.35 h Mpc~^ increases the parameter (Tq given by (p!9|), which strengthens the 
efficiency of the non-Gaussianity. In particular, the high frequency cutoff with the initial skewness 
hi = 3.0 has the curious behavior that the present correlation coefficient becomes larger than 
unity, which cannot be allowed because the Schwarz inequality states — 1 < r^. < 1. As for the 
non- vanishing initial kurtosis /12, we obtain the same present values ht and as obtained in the 
cutoff kc = 27r/8 h Mpc~^ case. However, the effective biasing induced by the non-Gaussianity is 
significant and we have the small value of the initial parameter 69 = 2.42. 



3.3.2. Quasi-linear scale 

On quasi-linear scales A; ^ 0.1 h Mpc^^, the next-to-leading order of non-linear gravity gives 
the important contribution to the growth of total matter power spectra (see Figure 1). 

We pay attention to the non-linear growth, which may be related to the scale-dependence of 
the biasing shown in Figure 2. From the weakly non-linear analysis Pmik) by several authors, we 
have recognized that the mode-mode coupling of long-wave modes contributes to the amplitude 



of the short-wavelength fluctuations (Makino, Sasaki & Suto 1992). This leads that the power 
spectrum evolves with the different growth rate, depending on the effective spectral index n^ff. 
For the typical spectrum like the CDM model, the one-loop correction causes the significant 
enhancement of the high-A; part of the power spectrum ( pain fc Bertschinger 199^ ). 

The weakly non-linear effect also appears in the power spectrum of galaxies Pg{k) and the 
cross correlation Px{k) (see eqs.|A-£]-[ A-12|1 ). On appropriately small scales, we could find the 
enhancement of perturbation 5g, whose growth rate depends on the initial parameters bQ,rQ and 
the initial non-Gaussianity, in addition to the effective index n^ff. Hence, the relative differences 
between the growth rate of 5g and 5m dynamically lead to the scale-dependence of the biasing (see 
Figure 2). 

To see the effect of non-linear growth explicitly, we write the parameters and as 



doop 



and Tfc 



41) 



loop doop 



and are the one-loop correction for the biasing 



parameter and the correlation coefficient. They are given by 



Aoap 



loop 



2 



P. 



^''\k) + P^''\k) 



P^'\k) + P^'\k) 



pf\k)+pr\k) 
p(")(fc) 



(13), 



Pm (A^) 
,(22) 



pr'{k)+p^'''\k) 

Pt\k) 



+ 



pis^\k) + p^^\k) 



p{ll) 



(k) 



The scale-dependence induced by gravity can be observed from the time evolution of b''^"^ and 

loop 



- 13 - 



In Figure 4, the snapshots of b^"^ and r^"^ are shown as the function of k for the different 
redshift z. The results at redshift z = are also displayed to emphasize the qualitative behavior. 
In any cases, the scale-dependence of the biasing cannot become significant within the validity 
of the perturbation theory, however, the scale-dependence obviously exists. The variation of the 
parameters })^°^ and r^""^ becomes larger on smaller scales. 

The characteristic behavior can be recognized in the two cases 60^0 > 1 and 60^0 < 1) which 
can be deduced from the one-loop spectrum (|A-9|) - (|A-12| ). For our qualitative understanding, we 
simply restrict the analysis on the /ii = /12 = cases. Figure 4a is the case 60^0 > li which has the 
initial condition bk = 3.0, = 0.8 at z = 5. Compared to the total mass distribution, the strong 
clustering of initial galaxies well-correlating with the total mass enhances the non-linear growth of 
5g. Thus, the one-loop correction of the biasing parameter and the correlation coefficient increases 
in time. On the other hand, the initial distribution 5g with the low correlation does not lead to 
the enhancement of clustering feature, but rather weakens the growth of fluctuations to have the 
sufficient correlation with the total mass. Figure 4b corresponds to this situation, in which we 
have chosen the condition as b^ = 3.0, = 0.2. Because the growth of perturbation 5ra becomes 
larger than that of 5g, the biasing parameter 6^^°°^ decreases. Interestingly, the one- loop correction 
rl^°^ initially decreases, although the evolution eventually strengthen the correlation. 



4. Summary and Discussion 

In this paper, we have focused on the gravitational evolution of the stochastic galaxy biasing 
deduced from the power spectrum. Using the perturbation theory, we develop the analysis taking 
into account the one-loop correction of the non-linear gravity. Assuming the non-local distribution 
/(Am.) for the total mass and the local initial condition ^(Ag) for the galaxy distribution, 
we examine the weakly non-linear power spectra and investigate the time evolution of biasing 
parameter and the correlation coefficient. The important conclusions are summarized as follows: 

• We found that the initial non-Gaussianity affects the galaxy biasing significantly. In the 
deterministic biasing, the time evolution of the biasing parameter with the non-vanishing 
initial skewness deviates from the linear theory prediction by Tegmark & Peebles (1998), 
even in the linear regime. As time goes on, the deviation cannot be neglected. However, the 
deviation can be reduced in the presence of the initial stochasticity. On the other hand, the 
correlation coefficient does not suffer from the initial non-Gaussianity. The time evolution of 
the correlation coefficient is well- approximated by the linear evolution. 

• The local initial condition of the galaxy distribution leads to the almost scale-independent 
evolution of the biasing on large scales. However, on quasi-linear scales, the one-loop 
correction of gravity becomes important and it can give the scale-dependence of the biasing. 
In the typical cases with 60 ''o > 1 and b^r^ < 1, we have the different prediction of the 
scale-dependence. For the low correlation b^r^ < 1, we have the small biasing parameter, 
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compared with the one on large scales. 



The prediction of the scale-dependent biasing just comes from the weakly non-linear results. 
We cannot say anything about the strong non-linear regime beyond the one-loop correction. 
However, the weakly non-linear analysis is suitable in investigating the transition from the linear 
to the non-linear regime. Therefore, our qualitative prediction will not change as long as the 
gravity plays a role for the scale-dependent biasing. In fact, the halo biasing property taking 
into account the merging process has recently shown that the biasing parameter decreases as the 
smoothing scale of density fluctuation becomes small. ( Mo fc White 199^ , Sheth fc Lemson 199£ , 



Jing 1998 ). Further, there is a suggestion that the halo-halo correlation can become anti-biasing 



on the relatively large scale A; ~ 0.15 h Mpc ^ ( Kravtsov &: Klypin 199S| ). This indicates that the 



gravitational clustering gives an important contribution. Therefore, the qualitative behavior can 
be partially understood from our weakly non-linear analysis. 

As a demonstration, we plot the power spectra Pm{k) and Pg{k) at z = 2.2 and 0.0 in 
Figure 5 {left panel). The parameters are given at z = 2.2, chosen as = 3.5, = 0.1 and 
/ii = /i2 = 0.0, i.e, we have the large scatter in the initial 6g-5m relation. On the scales k ^ 0.6 h 
Mpc~^, the perturbation theory breaks down at z = 0.0 that the one-loop correction of the total 
mass dominates the linear order power spectrum Pm^\ However, within the validity domain, 
the weakly non- linear spectrum Pg{k) decreases its power in the high frequency part and the 
anti-biasing property appears. This behavior can be illuminated by comparing with the N-body 
simulation in the right panel of Figure 5, where the spectrum of dark matter particles Pm{k) and 
halos Ph{k) are presented. The halo and the dark matter catalogue of the simulation are based 
on the standard CDM model ( Magira, Jing fc Suto 2000| ). The halos are identified using the 



Friend-of- Friends algorithm and selected with the mass threshold Mth = 4.47 x lO^^M©. Although 
the halo biasing at z = 2.2 is rather larger than the weakly non- linear biasing, the similarities 
at z = 0.0 is apparent. For the detailed comparison, further implementation of our analysis is 
needed, including the merging process of the halos. The volume exclusion of the halos could be 
essential for the anti-biasing on small scales. However, the similarities seen in Figure 5 implies 
that the gravitational evolution and the stochasticity of the galaxy biasing play important role for 
the present biasing features. We will proceed to improve our treatment quantitatively to explain 
the halo biasing. 

Throughout this paper, we have neglected the galaxy formation. For a realistic situation 
including this process, the cosmological simulation by Blanton, Cen & Tegmark (1999) shows 
that the hot gas in the overdense region prevents the star formation and the relation between the 
galaxies and the dark matter becomes a poor correlation. The realistic process will play a role for 
explaining the origin of stochasticity. Combining it with our analysis might yield the interesting 
results. 



Potentially, the effective biasing parameter due to the initial non-Gaussianity can cause a 
serious problem. The determination of the parameters /ii, /i2 and the correlation coefficient 
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from the observational data would become difficult if there exists the significant contribution of 
the non-Gaussianity to the galaxy clustering. The method using the higher order correlation 
1994, Paztahaga &: Frieman 1994 ) should be used carefully, which is partly mentioned by Taruya 



& Soda (1999). This might lead to a more practical issue, that is, the theoretical prediction of 
the high redshift objects such as the quasars. In this case, the various cosmological effects should 
be taken into account. The cosmological redshift distortion and the light-cone effect appeared 



in the statistical quantities are not negligible for the data analysis ( Matsubara &: Suto 1996 



Yamamoto &: Suto 1999 ). The situation will become more complicated in the case of the biasing 



being stochastic and non-linear. We shall clarify this issue in the next task. 

The author would like to thank Y.Suto for careful reading of the manuscript and stimulating 
discussions, and J. Soda for suggesting the idea of this work. We are also grateful to H.Magira and 
Y.P. Jing for providing us the data of the simulation catalog. 
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Appendix: Calculation of the one-loop power spectra 

In this appendix, we compute tiie one-loop correction of the power spectra. 

We start to obtain the second and the third order solutions of 6m and 6g. To get the solutions, 
it is convenient to write down the Fourier transform of the evolution equations (P), (0) and (^. 
The quantities 6m, 6g and 6 are expanded as 

_ , , f d k ikx c /I \ f\ { \ f ^ ^ ikx c /I 

6m,g[X,tj = J ^2^j3 ^ 6m,g[k,t), u[X,t) = J ^^^^3 ^ 6m,g[k,t), 

Then the evolution equations are rewritten with 

aMM + e(k.t) = - j ^v{k\ k - k')Uk - k', t)e(k', t), 

= - / TTTT^ [^(^'' k-k')+ V{k - k', k') - 2£{k', k - k')] 9{k - k' , t)e{k\ t), 

^ ^^^^ _ k%ik - fe', mk', t), 



where we define 



V{ki,k2) = 1 + \ \ C{ki,k2) = l 



Substituting the linear order solutions (p!o|) into the right hand side of the above equations, the 



second order solutions satisfying the initial conditions become (Pry 1984, Taruya, Koyama & Soda 



(A-1) 



19991) 

= J ^-^^ Snik - fci - k2) Am{ki)Am{k2) [D\t){Ti,2 - ^A,2) + lE{t)C,,2] 

5f){k,t) = 6(^Hk,t)-6(i\k,U) 

+ jfh^SD{k-ki-k2){Ag{ki)Ag{k2)-{27rf{Al)6D{ki)6D{k2)} 

_ fe, _ fc2)Pi,2 \Am{k^)Ag{k2) - A„(fci)A™(fe2)} , 

(A-2) 

where the expressions V(k\,k2) and C{ki,k2) are respectively suppressed as Vi^2 and £1,2- 
The solutions (|A-1| ) and ( [A-2[) contain the function -E(t) which satisfies E{ti) = 1. This is the 
inhomogeneous solution of the following equation: 

E + 2HE - -H'^nE = SH'^nD'^ + -D'^. 
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In Einstein-de Sitter universe, we have 



It is known that the $7 and A dependences of E/D"^ is extremely weak ( [Bernardeau 199^ ). 
Therefore, we develop the third order perturbations by replacing E{t) with (34/21)D^(t). Then 
we obtain the solutions 



7^7^1,2 (7'3,12 + 7^12,3) + l{Vl,2 - ^-Cl,2) (7^12,3 + 7^3,12 " '^J^12,3) 
-iA,2(P3,12 + 2Pl2,3) 



(A-3) 



k,t) 5^ij}{k,ti 



+ —h2 
24 



d'kid^k^dkl ^^^^ _k,-k2- fea) A,{k,)A,ik2)A,ik,) 



+ 



d?kidPk2dkl 

j—^ — ^ Soik -ki-k2- fea) 

7^3,12 {-^PW - 1) (^1,2 + 7^2,1 - ^^,2) A™(fci)A„(fc2)A™(fc3) 

+]^{D{t) - l)2pi,2 (A„(fci)A<,(fc2) - A„(fci)A„(fc2)) A„(fc3) 

+^hi {D{t) - 1) (A^(fci)A^(fc2) - {27rf{Al)6Diki)SD{k2)) A„(fc; 

+ lri2,3iD\t) - 1) (Vi,2 + V2,i - ^£1,2) Am{ki)Amik2) {Agiks) - A„(fc 

(A-4) 

where the quantities Vij^i, Vi^ij and Cij^i are respectively given by 

Vijj = V{ki + kj,ki), Vi^ij = 'P{ki,ki + kj), dj^i = C{ki + kj ,ki). 

Using the above results, the one-loop correction for the power spectra can be obtained by taking 
the ensemble average (H): 



Total mass distribution 



The one-loop power spectra Pm (k) and Pm^\k) are evaluated from the ensemble average 

= {27rf6D{k, + k2) P^^Hki), (A-5) 
2(<5«(fci)43)(fc2)) = {27rf6D{k, + k2) Pj^^\k,). (A-6) 



,(13) 
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The results are 



2 



_fq_ 

(2^ 



Pliq)Pl{\k - q\) a\q,k-q) 



p(^'Hk) = 2D' mk)J -^Pi{q) A{k,q,-q), 



(2vr)= 

which coincide with the weakly non-linear analysis of Makino, Sasaki & Suto (1992), Jain 
Bertschinger (1994) and Scoccimarro & Frieman (1996b). 

Galaxy distribution 



(A-7) 
(A-8) 



The power spectra P^^^\k) and P^^''^\k) are also given by the ensemble averages ( [A-5| ) and ( [A- 61 ) 
replacing the subscript (m) with (g). The resulting expressions are 

d^q 



:Pl{q)Pl{\k - q\ 



{D^ - If cx\q, k~q) + {D'-~ 1){D - l){b^r^ - 1) a{q, k - q)(3{q, k ~ q) 



+ {D - lYiboro - ly l{q, k-q) + {D- 1)'% - 26oro + 1) e{q, k - q) 
+ \i{D^ - l)(6or-o)' a(q, k - q) + hii{D - l)6gro(6o - ^o) /3(q, k - q) 

+ — hl bl 
18 ^ ° 



(A-9) 



P(^^\k) 



2 Pi{k) 



d^q 

(2^ 



(D-1 + boro) {(£>' - 1) A{k, q, -q) - {D - I) B{k, q, -q) 

-{D - if C{k, q, -q) + {D^ - l)(6oro - l)D{k, q, -q 
+ {D - if { (boro{D - 1) + bl) C{k, q, -q) + bl{rl - l)E{k, q, -q)} 

+ \\h boro{D - 1) + ^2 bl\ {{D - l)6oro + . (A-10) 



Cross correlation 

The one-loop power spectra for the cross correlation are computed from 

(5f (fcl)<5(^)(fc2)) = {27Tf5D{k,+k2)P^^'\ki), 

{6^J,Hk,)5^^Hk2) + 6l'\k^)6l^\k2)) = (2vr)3fe(fci + fc2)Pf^(^i). 

Then we have 
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{D^ - 1) k-q) + {D- l)(6oro - 1) /3(q, k - q) + -hi {h^roY 



(A-11) 



D{D^ - 1) + D^h^ro + D-1)) A{k, q, -q) - D{D - l)B{k, q, -q) 
+ (6oro - 1)d[{D - lfC{k, q, -q) + {D^ - l)D{k, q, -q)} 



+D boro |^/ii6oro(Z) - 1) + ^/is^o} 



(A-12) 



The bold-math quantities appeared in the expressions ( |A-7| )-( X-12D are defined as follows: 

4 

a{ki,k2) = 7^1,2 + 7'2,i - i^>Ci,2, 

/3(fci,fc2) = 7^1,2+^2,1, 
7(fcl,fc2) = 7^1,27^2,1, 

e(fci,fc2) = ^(7^1% + 7^1,1), 

A(fci,fc2,fc3) = ^ |7'3,12 (7^1,2 +7^2,1-^^,2) +7^12,3 (7^1,2 +7'2,1-^A,2) I 

+ ^ (7'l2,3 + 7'3,12 - 2Ci2,3) (ri,2 + 7^2,1 " ^Ci,2 
B{ki,k2,k3) = 7'3,12 (7^1,2 + 7^2,1 - ^/:i,2) , 
C{ki,k2,k3) = ^7^3,12 (7^1,2+ 7^2,1 ), 
D{ki,k2,ks) = ^7^12,3 (7'l,2 + 7^2,1 - ^A,2 
E{ki,k2,k^) = ^7^l,27'3,12. 

To develop the further evaluation of the one-loop power spectra, we write the integrals in spherical 
coordinates q, 1!) and (p, where we choose the polar axis aligned with the wave vector k. Then 
we introduce the new variables x = q/k and u = cos-d. The integral over the azimuthal angle (p 
simply yields the multiplicative factor 2tt to each one-loop spectrum. As for the integration of u, 
the P^^^^(fe) part needs the explicit expression of the linear spectrum Pi{k), while we can integrate 
the part without knowing the function Pi{k). Thus, using the explicit form Pi{k) given by 

(pOD, we numerically evaluate the integral over x and u for the P^^'^^k) part. For the P'^^^^k) 
part, the one-dimensional integration over x is only required. 
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Table 1. The list of model parameters used in Figure 3 and the initial conditions bk, Vk given at 
z = h for the fixed mode k = 27r/100 h Mpc~^. Together with the linear results 6^^^ and r^^\ the 
present parameters and are also given by extrapolating the weakly non-linear analysis. 



cutoff length 




model parameters 




z ■ 


= 5 


z 


= 


2iT/kc [h-^ Mpc] 


bo 




hi 


h2 


h 


rk 


bk {b(^^) 


Tk (r-(')) 


8 


3.0 


0.8 


0.0 


0.0 


3.00 


0.80 


1.44(1.27) 


0.99(0.97) 


8 


2.93 


0.82 


3.0 


0.0 


3.00 


0.80 


1.69(1.27) 


1.00(0.97) 


8 


2.93 


0.82 


-3.0 


0.0 


3.00 


0.80 


1.13(1.27) 


1.00(0.97) 


8 


2.64 


0.79 


0.0 


3.0 


3.00 


0.80 


1.40(1.27) 


0.99(0.97) 


8 


3.0 


0.20 


0.0 


0.0 


3.00 


0.20 


1.03(1.05) 


0.85(0.89) 


8 


2.93 


0.205 


3.0 


0.0 


3.00 


0.20 


1.08(1.05) 


0.85(0.89) 


8 


2.93 


0.205 


-3.0 


0.0 


3.00 


0.20 


0.97(1.05) 


0.91(0.89) 


8 


2.6i 


0.199 


0.0 


3.0 


3.00 


0.20 


1.02(1.05) 


0.88(0.89) 


8 


0.50 


0.80 


3.0 


0.0 


3.00 


0.80 


0.85(0.91) 


1.00(1.00) 


8 


0.50 


0.80 


0.0 


3.0 


3.00 


0.80 


0.85(0.91) 


1.00(1.00) 


4.35 


2.91 


0.824 


3.0 


0.0 


3.00 


0.80 


1.97(1.27) 


1.05(0.97) 


4.35 


2.42 


0.782 


0.0 


3.0 


3.00 


0.80 


1.41(1.27) 


0.99(0.97) 
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b^^ = 3.0, r^ = 0.8 at z = 5 




1 1 1 1 I I I I I I 

0.01 0.1 

k [h Mpc-i] 



Fig. 1. — The Evolved results of power spectra at z = 0. The thin-dashed line shows the linear 
power spectrum of the total mass Pm^\k) evolved from z = 5. The initial spectrum is given by the 
CDM spectrum (20) with the appropriate model parameters (see text). The thin-solid, thick-solid 
and thick-dashed lines are the weakly non- linear power spectrum of the total mass, galaxies and 
the cross correlation between them, respectively. The initial conditions at z = 5 are specified as 
bo = 3.0, ro = 0.8 and hi = h2 = 0. 
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Fig. 2. — The snapshots of the evolved biasing parameter (upper panel) and the correlation 
coefHcient (/oxyer panel) are depicted as the function of the Fourier mode k. The initial condition 
are the same as in Figure 1. Within the interval 0.001 < k < 0.1 h Mpc~"^, the biasing parameter 
and the correlation coefficient becomes almost independent of scales. 
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bj^ = 3.0, at z = 5 




I I I I I I I I I I I 

5 4 3 2 10 

redshift z 



Fig. 3a. — Time evolution of the biasing parameter bk {upper panel) and the correlation coefficient 
rk {lower panel) as the function of the redshift z, fixing the Fourier mode k = 27r/100 h Mpc~^. By 
adjusting the model parameters bo and ro, the initial conditions at z = 5 are specified as b^ = 3.0, 
Tfc = 0.8, i.e, almost deterministic biasing case. The initial non-Gaussianity are respectively chosen 
as /ii = /i2 = 0.0{solid line), hi = 3.0, /12 = 0.0{long- dashed line), hi = — 3.0, /12 = 0.0 {short- dashed 
line) and hi = 0.0, /12 = 3.0{dotted line). The results are extrapolated to the present time z = 0. 
For comparison, we also plot the linear result given by ( |l6| ) with the same initial condition bo = 3.0, 
ro = 0.8 {thin-solid line). 
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bj^ = 3.0, at z = 5 
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Fig. 3b. — The same figure as in Figure 3a, but with the different initial conditions 
rfc = 0.2 {stochastic biasing case). 
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bj^=3.0, r^=O.Q at z = 5 
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Fig. 4a. — The one loop contribution of the biasing parameter bj^"^ and the correlation coefficient 
r^""^ as the function of the Fourier mode k. The initial conditions given at 2; = 5 are specified 
as 6jt = 3.0 and = 0.8. For simplicity, the initial non-Gaussianity are set by h\ = h2 = 0. 
Because of the cutoff kc, the initial parameters and r^°^ slightly lose their power at high 
frequency part. However, it is obvious that the spatial dependence of the parameters 6^'°"^^ and 
^(loop) (jyj^g^j^jpg^^jy appears. On smaller scales, the spatial variation becomes larger. The typical 
behaviors are classified as 6oro > 1 (Figure 4a) and hoVQ < 1 (Figure 4b). 
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bj^=3.0, r^=0.2 at z = 5 
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Fig. 4b. — The same figure as depicted in Figure 4a, but with the initial conditions bk = 3.0 and 
r/j = 0.2, corresponding to the 6oro < 1 case. 



-29- 



n 
O 



CO 

I 



CO 

O 



CO 

I 



Weakly non-linear 



I I I 1 1 1 I I I n= 



^ 1000 h- 



100 r 



^ 1000 



100 - 



10 



10 - I I I III 



z = 2.2- 




I I I I I- 
z = 0.0 = 




I I I I I L 



Simulation 

— I — I — r~r 

.z = 2.2 



I- 1 1 1 1 1 1 1 1 — I — I I i-i 



- I 1 1 III 



I I I I 
z = 0.0^ 




Ph(k) 



I I I I I L 



0.1 

k [h Mpc-i] 



0.1 

k [h Mpc-i] 



Fig. 5. — The anti-biasing property in the weakly non-hnear analysis (left panel) and the N-body 
simulation (right panel). In left panel, the weakly non-linear power spectra Pm{k) and Pg{k) are 
depicted under the initial conditions = 3.5, = 0.1 and hi = 1x2 = 0.0 given at 2; = 2.2. In 
right panel, using the simulation catalog of Magira, Jing k. Suto (2000), the numerical results of 
dark matter particles and the halos are presented in the case of Standard CDM model. The halos 
are identified using the FOF algorithm and selected with the mass threshold Mth = 4.47 x W^'^Mq. 
Qualitatively, the weakly non-linear result can recover the anti-biasing feature seen in the simulation 
if the large scatter in the 5g-5m relation is present initially. 



